Mesoscopic phenomena in Bose-Einstein systems: 
Persistent currents, population oscillations and quantal phases 
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Mesoscopic phenomena — including population oscillations and persistent currents driven by quan- 
tal phases — are explored theoretically in the context of multiply-connected Bose-Einstein systems 
composed of trapped alkali-metal gas atoms. These atomic phenomena are bosonic analogues of 
electronic persistent currents in normal metals and Little-Parks oscillations in superconductors. 
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Introduction: The purpose of this Letter is to consider 
multiply-connected many-particle systems obeying Bose- 
Einstein statistics and, in particular, to address the sen- 
sitivity of such systems to quantal phases. If the bosonic 
constituents (e.g., alkali- metal gas atoms) are electrically 
neutral, electromagnetism, in the form of the Aharonov- 
Bohm (AB) phase, cannot provide a source of quantal 
phase. Then, in seeking sensitivity to quantal phases 
we are led to consider the spin degree of freedom of 
the bosons, and the consequent possibility of quantal 
phases of geometric origin pj . As we shall see, geometric 
quantal phases can readily affect the energy levels, and 
hence populations, of the single-particle quantum states, 
and lead to persistent equilibrium currents in multiply- 
connected systems, thus providing a striking example of 
quantal mesoscopic phenomena in the setting of bosonic 
systems. Such phenomena are bosonic analogues of phe- 
nomena well known in the context of the mesoscopic 
physics of normal-state electronic systems (such as per- 
sistent equilibrium currents and conductance oscillations 
in conducting rings) which arise due to AB [^j or geomet- 
ric H quantal phases. They are also bosonic analogues 
of the flux-sensitivity of the superconducting transition 
temperature of a thin superconducting ring, known as 
Little-Parks oscillations Q. (These quantum interfer- 
ence phenomena are mesoscopic, in the sense that they 
vanish in the limit of large system-size.) 

It is worth mentioning that there is a sense in which 
bosonic settings are preferable to electronic settings, if 
one wishes to observe implications of quantal phases in 
many-particle physics: in the fermionic case, the Pauli 
exclusion principle forces the occupation of many single- 
particle states, and there are strong cancellations be- 
tween the effects of quantal phases on these states. By 
contrast, Bose-Einstein statistics promote the signifi- 
cance of the single-particle ground state. In this sense 
then, bosonic systems tend to amplify mesoscopic effects, 
at least in comparison with fermionic systems. 

One scheme for introducing a geometric quantal phase 
is to have the bosons move through regions of space in 
which there is a spatially varying magnetic field to which 
the spins of the bosons are Zeeman-coupled. Then, as 



discussed in Ref. in the context of magnetic traps, 
the inhomogeneous magnetic field (if sufficiently strong) 
leads to a geometric vector potential A (and a corre- 
sponding geometric flux $) which influences the orbital 
motion of the bosons, and does so in much the same way 
as the electromagnetic (AB) vector potential (and flux) 
influences the motion of electrically charged particles. 




FIG. 1. Toroidal sample of non-circular cross-section. 

References || considered conventional (i.e., not AB- 
like) consequences of the geometric vector potential, i.e., 
effects associated with nonzero values of the geometric 
field-strength Jl = V x A (i.e., the vorticity). However — 
and this is the main point of the present Letter — there 
are striking quantal AB-like consequences of the geomet- 
ric vector potential itself (rather than the field-strength) , 
especially in multiply-connected configurations, and even 
when fi = in the sample. These consequences include 
an oscillatory dependence on the geometric flux <!> of the 
energies of all single-particle levels and, thus, the equi- 
librium populations of these levels and, more generally, 
all equilibrium quantities. Moreover, these single-particle 
energy-level oscillations lead directly to the existence of 
equilibrium currents |(| that flow around the trap (at 
generic values of <I>). Therefore, despite the electrical neu- 
trality of the atoms in BEC systems, the geometric phase 
allows one to realize counterparts of the well-established 
collection of electromagnetic AB-phase sensitive phenom- 
ena, as well as new, bosonic, phenomena such as large os- 
cillations in the populations of the various single-particle 
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levels. 

In order to realize geometric-phase-driven oscillations 
one needs to vary the inhomogeneous magnetic field |t]] 
that the bosons inhabit. In magnetic traps, the abil- 
ity to make such variations is limited, as variations alter 
the structure of the trap (i.e., the shape of the system). 
The recent achievement of confinement via purely opti- 
cal traps || liberates the magnetic field from its dual role 
of confining the atoms and causing the geometric phase, 
and thus enlarges the scope for exploring the effects of ge- 
ometric phases. We shall touch upon the issue of making 
multiply-connected traps at the end of this Letter. 
Model system: We consider a system of many identi- 
cal charge-neutral noninterating bosonic atoms of mass 
m, confined to a multiply-connected trap. For the sake 
of simplicity, we envisage the trap as having the follow- 
ing features: (i) it is toroidal and axisymmetric; (ii) it 
is sufficiently narrow in the radial direction that, under 
operating conditions, only the state with the lowest ra- 
dial quantum number is occupied (i.e., the radial energy 
scale TiQ, r is large) ; (iii) in the axial direction there is con- 
finement by an oscillator potential; and (iv) if the trap is 
optical then only mild conditions need be obeyed if an ap- 
plied inhomogeneous magnetic field is not to change the 
trapping potential appreciably, (v) If the trap is mag- 
netic then the application of an additional homogeneous 
field will change $ but not the confining potential. In 
the absence of <£, the spectrum of single-particle energy 
eigenvalues £® n is given by £ ® n = hfl z n + Tif^ £ 2 , where 
h£l z is the axial oscillator energy scale, S.f^ is the az- 
imuthal energy scale, we have omitted all zero-point en- 
ergy contributions, and we do not allow for radial excita- 
tion. The quantum numbers £ and n range, respectively, 
over all integers and all non-negative integers. 

We now suppose that the toroidal system of trapped 
atoms is subjected to a magnetic field having the prop- 
erty that the orientation of the field varies across the 
trap ||. Under conditions of adiabaticity (i.e., f2^ much 
smaller than the Larmor frequency of the spins) for the 
dynamics of the spins of the atoms, the dominant effect 
of the inhomogeneous magnetic field on the energy spec- 
trum is to introduce a spin-dependent flux $, so that the 
spectrum (for atoms with spin-projection lying parallel 
to the magnetic field) becomes 

£ t , n = hn z n + nn^(£-<t>) 2 . (i) 

(Atoms with spin-projection lying antiparallel to the 
magnetic field direction are not trapped.) The origin of 
this flux is the Berry phase associated with the spin dy- 
namics (see, e.g., Refs. Our aim is to compute the 
number of particles in the single-particle ground state as 
a function of the temperature T, the geometric flux 
and the mean total number of particles N, and to do 
so via the grand canonical ensemble Jig ], To this end, 
we first compute the total number of particles N, as a 
function of T, $ and the chemical potential \i: 

]T Nt, n , (2a) 

£— — 00 * — ^n— 



N t , n = { e (^»-f)/W-i}- 1 . (2b) 
Next, we decompose this sum: N = N + iV, where 

To compute N, we use a variant of the contour integra- 
tion technique described, e.g., in Refs. pj] ], which yields 

~ 7T 

N = {cot 7T ($ — <p) — COt 7T ($ + tp)} 

2tp 

+0(exp(-(2 7 r) 3 / 2 /^ /2 )), (4) 

where we have introduced the reduced frequency = 
Ml^/ksT, and the reduced chemical potential jl = 
and, for convenience, ip denotes (/t/^) 1/2 • [The 
form for N arrived at by this technique is much more 
rapidly convergent than the original form, Eq. (|2a|)1. Not 
surprisingly, however, even after omitting exponentially 
small terms, the transcendental equation for jl(T, $, N), 
Eq. (||), cannot, in general, be solved explicitly. With- 
out loss of generality, let us assume that < $ < 1/2. 
(Results for other values of $ can be obtained via the 
symmetries of reflection, $ — > — $, and translation, 
<f> — > $ + 1.) Then further simplification is achievable 
in three cases: (i) $ -C 1/2, (ii) $ > (M^)- 1 / 2 , and 
(iii) |cot27r$| 3> Nfl^/ir. In case (i) we expand the 
cotangents on the r.h.s of Eq. (||) in Laurent series, re- 
taining two terms in each series, and solve the resulting 
equation for jl, thus arriving at 

fl^n^ 2 -{N- (tt/3^)}- 1 . (5) 

In case (ii) we instead expand the prefactor and the ar- 
guments of the cotangents in Eq. (|j) to linear order in 
the deviation of jl from the (dimensionless) single-particle 
ground-state energy Cl^Q 2 . Thus, we arrive at 

jl « f2 $ 2 - 2tt~ 1 <I> cot -1 ( cot 2tt$ + Tr^O^iV 

+ v /cot 2 7r$ + ($^7V/7r) 2 + l V (6) 

In case (iii), which corresponds to <E> close to 1/2, we have 
jl w fi $ 2 - 2tt" 1 $ cot" 1 ($^A>/vr) . (7) 

In making our expansions of Eq. (Q) for N we restrict 
ourselves to the regime of quasi-BEC, i.e., we consider 
values of [i only slightly smaller than the single-particle 
ground-state energy £"0,0 (i.e., fiVL^ 2 ). 
Physical consequences of the geometric flux: To use these 
results for jl(T, N) to determine desired physical quan- 
tities, such as the populations N^ n of the single-particle 
states as functions of the variables (T, N), we first note 
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that in the regime of BEC we need only retain the dif- 
ference between fi and £q.o in N, but may omit 




Geometric flux 



FIG. 2. Oscillatory dependence on the geometric flux of the 
populations of the lowest three single-particle energy levels at 
a higher (dashed line) and a lower (full line) temperature. 

it from N. Next, we observe that N = N + N, so that 
by knowing ji(T,$,N) we know fl(T,$,N). This we 
use to eliminate fi from the Bose functions that deter- 
mine the populations iVf )n . To illustrate the oscillatory 
behavior of single-particle state populations with $ we 
show, in Fig. ||, the populations Ni l7n of the three lowest- 
lying states (i.e., Nq.q, 2Vi,o and iV 2 ,o, as |$| < 1/2) as 
functions of <& |lSq| . We have chosen for this illustra- 
tion a system of N = 10 5 atoms of 87 Rb, sample ra- 
dius R — 1 pan, axial trap frequency Q z = 7, 500 Hz, 
temperature T = 5^K (for which 800^ = 400* = 1) 
or T = 10/iK (for which 4OOfi = 20n z = 1), where 
Cl z =hQ z /k B T. Note that at $=0 we have N^ n =N-t n , 
and at $=1/2 we have Ne tn =Ng+i in . The latter case 
illustrates the more general point that at half-integral 
values of $ the lowest single-particle energy level is de- 
generate for the case of perfectly azimuthally symmetric 
traps. If, however, the azimuthal symmetry is absent 
then the level crossing is avoided, and the single-particle 
ground state is separated from the excited states by an 
energy gap at all values of $. 

The population oscillations are mesoscopic, in the 
sense that they vanish in the thermodynamic limit [ p^[ . 
Indeed, the number of atoms in traps, although typi- 
cally large, is not on the order of Avogadro's number 
and, therefore, the systems are even further from the 
thermodynamic limit than conventional macroscopic and 
even mesoscopic samples. Nevertheless, although BEC is 
not, strictly speaking, marked by a sharp thermodynamic 
phase transition, atomic condensates acquire features of 
the thermodynamic limit already at N = 10 4 (see, e.g., 
Refs. [ fjr|]l5"| ) . Despite this, our calculations reveal oscil- 
latory phenomena at N = 10 5 whenever the fraction of 
atoms in the ground state is appreciable. Thus, one has 
the capability of observing, simultaneously, both macro- 
scopic and mesoscopic phenomena. 



The ^-dependence of the single-particle energy lev- 
els also leads to the phenomenon of equilibrium persis- 
tent currents, i.e., dissipationless particle-currents that 
flow around the trap. (Such equilibrium currents should, 
of course, be distinguished from the noncquilibrium 
metastable currents that can arise in multiply-connected 
samples; see, e.g., Refs. jj.) The simple formula, 
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FIG. 3. Oscillatory dependence of the equilibrium persis- 
tent particle-current (per atom) on the geometric flux at a 
higher (dashed line) and a lower (full line) temperature. 

for the persistent particle-current / follows from the gen- 
eral relation 5F/5A = — j, where F is the free energy, 
A is a gauge potential (such as the geometric vector po- 
tential) and j is the conjugate current-density. In Fig. || 
we show the dependence of I/N on <E> at the conditions 
and temperatures specified in the previous paragraph. 
As T is increased from zero the saw-tooth form of the <&- 
dependence is smoothed to a more sinusoidal form, but 
the zeros at integral and half-integral values of the flux 
are preserved. As for the characteristic scale of the am- 
plitude of /, it is on the order of Ntt^, at least when the 
ground state contains most of the particles. 

0.42i . . 1 
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FIG. 4. Oscillatory dependence of the specific heat per 
atom on the geometric flux. 

As it is Fermi rather than Bose systems that have tra- 
ditionally provided settings for mesoscopic physics, we 
pause to compare the characteristic magnitudes of per- 
sistent currents in Fermi and Bose systems. In the special 
case of single-channel systems the characteristic magni- 
tudes are similar: the many bosons in the (low-velocity) 
ground state contributing roughly as much as the sin- 
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gle (high-velocity) fermion at the Fermi level (contribu- 
tions from fermions below the Fermi level essentially can- 
celling one another). In the more general case, however, 
in which there are many channels, the greater the extent 
of transverse excitation, the smaller the contribution to 
the persistent fermion current (owing to the reduced ki- 
netic energy at the Fermi level). By contrast, for bosons 
the particle occupations are, of course, not spread over 
the many current-reduced channels, and instead are con- 
centrated on the optimal channel. Thus, for bosonic sys- 
tems such mesoscopic effects are amplified, relative to the 
Fermi case. 

As a third consequence of the geometric flux, we con- 
sider the oscillatory behavior of the (dimensionless) spe- 
cific heat (per particle) C = dE(T,<f>, N)/dNk B T. (We 
recognize that C may be difficult to measure.) For C, 
a more pronounced ^-dependence arises in the case of 
traps that are weakly confining in the axial direction (i.e., 
Q z <C fi^). In Fig. f| we show C($) for the case of 23 Na 
(the reduced mass of which also enhances the sensitivity 
to $ compared with 87 Rb) at R = 0.5 ^m, T = 40 nK 
and Q z = 100 Hz. Upon closer inspection, the apparent 
cusp at the level crossing, which is a remnant of the true 
singularity at T = 0, is seen to be rounded. 
Experimental issues and concluding remarks: Having de- 
scribed several consequences of the geometric flux, we 
now briefly discuss some issues concerning the possibility 
of the experimental realization of these consequences. We 
see three pivotal matters: (i) how to construct a toroidal 
sample; (ii) how to subject the sample to a suitably in- 
homogeneous field; and (iii) how to detect population os- 
cillations and persistent currents. As for (i), it should be 
feasible to construct toroidal samples based on magnetic 
traps by using a blue-detuned laser to repel atoms from 
the trap center. We hope that in purely optical traps j8[ a 
similar method, combining red- and blue-detuned lasers, 
could be employed to make a toroidal sample. As for (ii), 
magnetic trap technology itself is suitable for creating 
the necessary textured magnetic fields. (As magnetic 
traps are currently used as a first stage in the loading 
of optical traps, creating textured fields should not im- 
pose a large additional experimental burden.) As for (iii), 
the phonon-imaging technique [|l6| discussed with regard 
to metastable currents in the last paragraph of the sec- 
ond of Refs. is less difficult in the present setting of 
equilibrium persistent currents, owing to the far larger 
magnitude of the latter. Other experimental possibil- 
ities may include (nondestructive) light-scattering and 
(destructive) timc-of-flight techniques. 

We note that interactions between the bosons are not 
expected to alter qualitatively the oscillatory effects that 
we have discussed. This is because, in the present con- 
text, k^T 3> /i, for which one expects interactions to have 
only perturbative consequences . We conclude by re- 
marking that if it is possible to realize superconductivity 
via the BEC of pre- formed bosons then oscillations in the 
level-populations, persistent currents, and perhaps spe- 
cific heat, similar to those described in the present Letter, 



may be observable as consequences of an AB flux. 
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